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Abstract 

Background: Extensive work has shown that vectors almost never feed at random. Often, a subset of individual 
hosts and host species are fed on much more frequently than expected from their abundance and this can amplify 
pathogen transmission. However, the drivers of variation in contact patterns between vectors and their hosts are 
not well understood, even in relatively well-studied systems such as West Nile virus (WNV). 

Methods: We compared roosting height and roost aggregation size of seven avian host species of WNV with 
patterns of host-seeking mosquito {Culex pipiens) abundance at communal and non-communal roost sites. 

Results: First, host-seeking mosquito abundance increased with height and paralleled increased mosquito feeding 
preferences on species roosting higher in the tree canopy. Second, there were several hundred-fold fewer 
mosquitoes per bird trapped at American robin {Turdus migratorius) communal roosts compared to non-communal 
roost sites, which could reduce transmission from and to this key amplifying host species. Third, seasonal changes 
in communal roost formation may partly explain observed seasonal changes in mosquito feeding patterns, 
including a decrease in feeding on communal roosting robins. 

Conclusions: These results illustrate how variation in habitat use by hosts and vectors and social aggregation by 
hosts influence vector-host interactions and link the behavioral ecology of birds and the transmission of vector-borne 
diseases to humans. 
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Background 

The behavioral ecology of hosts can play a significant 
role in determining pathogen transmission dynamics. 
For diseases where transmission is direct, host foraging 
ecology, habitat preference, and social interactions, in- 
cluding mating strategy, can influence the probability of 
contact with an infected host or environment and create 
hotspots for transmission [1-4]. However, the impact of 
animal behavior on the transmission of vector-borne 
pathogens is less clear [5,6]. Most studies exploring the 
influence of host behavior on vector contact rates have 
examined how human behavior influences the transmis- 
sion of pathogens causing malaria, dengue fever, and 
Chagas disease [7-9]. Understanding and controlling 
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infectious diseases that infect wild animals (many of 
which infect multiple species, including humans) re- 
quires examining the behaviors that influence contact 
between vectors and hosts [10-13]. 

An ideal free distribution model developed for host- 
seeking vectors predicts that vectors will feed on hosts 
to maximize their feeding success, and the degree of het- 
erogeneity among and within species will vary with host 
and vector densities, and host defensive behavior [14]. 
However, the importance of these factors, as well as host 
habitat use and aggregation in determining the extent to 
which vectors successfully move across the landscape to 
maximize feeding success is unknown. 

Although animals are well known to form communal 
groups to reduce the risk of predation [15], it is less 
clear whether aggregations form in response to biting ar- 
thropods and the pathogens these arthropods carry [16]. 
For example, although herd formation in caribou re- 
duced biting for individuals in the center of a herd, the 
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temporal pattern of formation of groups suggested ag- 
gregation was more for thermoregulatory reasons than 
in response to biting flies [17,18]. Similarly, while para- 
sitism has been shown to increase with sociality in some 
birds [19,20], in other taxa infection intensity of endo- 
and ectoparasites often decreases with increasing host 
group size [21]. At least as important as the influence of 
biting arthropods on the formation of communal groups 
is the reciprocal impact of host habitat use and group 
formation on the transmission of vector-borne diseases 
such as West Nile virus (WNV). 

West Nile virus was first detected in North America 
in 1999, and is primarily maintained in an enzootic 
cycle between avian hosts and mosquito vectors [11]. 
Nocturnal-feeding mosquitoes of the genus Culex, spe- 
cifically Cx. pipiens, Cx. restuans, and Cx. tarsalis are 
the most important enzootic and bridge vectors of 
WNV in urban, residential, and agricultural areas in 
North America [6,22-25]. Although over 200 species of 
birds have been found to be infected with WNV, 
mounting evidence from the mid- Atlantic, Northeast 
and Midwestern USA, and Colorado suggest that a 
single species, American robin (Turdus migratorius), is 
likely responsible for infecting the majority of WNV- 
infectious mosquitoes [26-31]. However, the mecha- 
nisms leading to preferential mosquito feeding on this 
species, and avoidance of others are unknown. 

One potentially important and largely unexplored as- 
pect of WNV transmission ecology is the behavior of 
hosts, including their nocturnal roosting patterns, which 
determine their location during the feeding periods of 
the dominant WNV vectors. Although host defensive 
behavior [32-34] and traits such as age [12,35] have been 
shown to influence the success of feeding mosquitoes, 
the influence of roosting behavior has rarely been ex- 
plored. Two aspects of nocturnal roosting behavior are 
likely to influence contact rates with foraging mosqui- 
toes. First, host-seeking Culex mosquito abundance may 
be higher in the forest canopy [36-39] so birds roosting 
higher may be exposed to more mosquito bites. Second, 
the formation of communal roosts could influence vector- 
host contact rates in several ways. Multiple birds produce 
more C0 2 than an individual bird, which could increase 
the number of mosquitoes attracted to the group. However, 
the increase in host seeking mosquitoes may not be propor- 
tional to group size, in which case communal roosting 
would decrease the average number of mosquitoes feeding 
on each bird. Previous research at communal roosts of 
American robins in Connecticut found fewer numbers of 
Culex mosquitoes were trapped at roosts versus non-roost 
sites, but WNV prevalence in mosquitoes was higher at 
communal roosts [40], whereas in Illinois there were no 
significant associations between distance to an American 
robin communal roost and WNV-infection in mosquitoes 



[41]. However, neither of these studies determined the 
number of host-seeking mosquitoes per bird at communal 
roosts and non-roost sites. The per capita number of 
vectors per host (i.e. vector-host ratio) is an important 
factor for vector-borne transmission and disease ex- 
posure [14,42,43]. 

We studied roosting height and roost aggregation size 
of seven bird species, and simultaneously quantified 
mosquito host-seeking abundance at the same sites. The 
seven species include the most important amplification 
host for WNV across a wide range of North America, 
the American robin [26-30,44], and the most abundant 
bird species in many urbanized and agricultural areas of 
North America, house sparrows and European starlings. 
We hypothesized that avian species roosting higher in 
the canopy and in smaller groups would have higher per 
capita contact rates with host-seeking mosquitoes and 
might explain preferential feeding by mosquitoes on those 
species [27,28,44]. We further hypothesized that mosquito 
abundance would be larger at communal roosts of birds, 
but not proportionate to roost size, such that the per capita 
number of mosquitoes would be smaller. In addition, be- 
cause large groups of birds emitting C0 2 could diminish 
the trapping efficiency of mosquito traps (which also use 
C0 2 as the primary attractant), resulting in biased estimates 
of mosquito abundance at communal roosts, we performed 
an experiment with artificial sources of C0 2 to estimate the 
magnitude of this effect. 

Methods 

Our study was conducted at three sites in Maryland and 
Washington, DC [25] including an urban area in Foggy 
Bottom, DC (38° 54' 07.33" N, 77° 3' 18.27" W), a resi- 
dential area in Takoma Park, MD (38° 58' 23.55" N, 77° 0' 
18.80" W), and a predominantly forested area in Rockville, 
MD (39° 6' 44.88" N, 77° 6' 23.80" W). These field sites 
were part of an ongoing multi-year study on WNV trans- 
mission in the greater Washington, DC area where data on 
mosquito feeding patterns, host abundance, and feeding 
preferences have been previously collected [28,29]. 
Data on avian nocturnal roosting behavior and host- 
seeking mosquito abundance were collected from July 
through September 2010, when WNV infection prevalence 
peaks in mosquitoes and when most human infections 
occur in the eastern USA [29,45]. 

Roosting behavior 

We attached transmitters (Lotek Wireless Inc. Newmarket, 
Ontario) weighing less than 5% of an individuals body mass 
[46] to seven bird species: American robin (transmitter 
weight: 2.8 g, N = 89 individual birds), Northern cardinal 
{Cardinalis cardinalis, 1.5 g, N = 60), tufted titmouse 
(Baeolophus bicolor, 0.7 g, N = 5), gray catbird {Dumetella 
carolinensis, 1.4 g, N = 15), mourning dove {Zenaida 
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macroura, 2.8 g, N = 8), European starling (Sturnus vul- 
garis, 2.8 g, N = 7), and house sparrow (0.7 g, N = 29). All 
birds except one American robin and six of the European 
starlings were juvenile (hatch-year) individuals. We focused 
on juvenile birds for the roosting studies because they ap- 
pear to play a key role in the transmission of WNV [47]. 
Birds were caught by removing nestlings from nests 
found by observing parental behavior (American robins 
and Northern cardinals), or were captured using mist- 
netting (all other species). Birds were then banded with 
an aluminum USFWS and unique color band combin- 
ation, radio-tagged, and released. All research with birds 
was approved by animal use and care protocol Kilpmlll2 
from the University of California Santa Cruz s Animal Care 
and Use Committee. 

We attached transmitters using a figure eight shaped 
two-loop harness design that slid over the legs of the 
bird and allowed the transmitter to rest on the syn- 
sacrum [48]. Harness loops were made from clothing 
elastic, which has been shown to work as a durable but 
less-constrictive material for short-term attachments of 
transmitters to birds [49]. Harness size was adjusted to 
individual body size and fixed with super glue to prevent 
unraveling. We observed post-banding movements of 
tagged individuals to insure attached transmitters were 
fitted properly and that nestlings continued to be fed by 
parents. 

We used Biotracker scanning receivers in combination 
and three-element folding yagi antennas (Lotek Wireless 
Inc. Newmarket, Ontario) to locate birds fitted with trans- 
mitters. To record nightly roosting location of tagged birds, 
bearings were taken simultaneously from two or more 
points for triangulation. Once the approximate location of 
a bird was discovered we used ToughCam thermal imaging 
cameras (Infrared Cameras Inc. Beaumont, TX, USA) in 
tandem with radio receivers to determine the roosting loca- 
tion and height of individuals [50]. We estimated the height 
of the bird using a clinometer, a 50 m tape, and the birds 
location. 

Over the course of the study we found several tagged 
individuals roosting in large communal aggregations 
while others remained solitary. We estimated communal 
roost aggregation size by counting the number of birds 
arriving from the four cardinal directions at dusk and by 
counting all birds within the contiguous vegetation area 
surrounding the roost with thermal imaging cameras. 
Roost size estimation methods were used concurrently and 
validated against one another to ensure both methods pro- 
duced similar results. The two American robin communal 
roost sites were 7,128 and 2,380 m 3 in size and at peak oc- 
cupation contained roughly 2.5 and 4.1 American robins 
per m 3 respectively. When a bird was found roosting soli- 
tarily, not in association with a large communal aggrega- 
tion, we searched within a -10 m radius of the bird using a 



thermal imaging camera to determine if additional hosts 
were roosting in close proximity. Previous research has 
shown that 10 m is the approximate distance that foraging 
mosquitoes respond to hosts [51]. 

Host-seeking mosquito abundance 

We measured the abundance of host-seeking Cx, pipiens 
mosquitoes, at four sites (two non-communal and two 
communal roost sites) in Takoma Park, MD using two 
to six CDC light traps hung at two heights: 1.5 and 12 
meters, and baited with C0 2 (-300 g dry ice per trap). 
We identified mosquitoes using morphological characteris- 
tics [52] and identified a subset of Culex mosquitoes using 
a molecular assay to determine the fraction of Culex mos- 
quitoes that were Cx. pipiens [53]. We were unable to test 
mosquitoes trapped at these sites for WNV because sam- 
ples were lost in shipment to the testing lab. 

Communal roost sites consisted of large American robin 
aggregations discovered while radio-tracking robins. Mos- 
quito traps at communal roost sites were placed within and 
on the edge of roosts. Non-communal roost sites were the 
location where birds were initially banded. All four sites 
(two communal and two non-communal) were within a 
circular area with an approximate radius of 10 km and 
consisted of relatively similar residential areas with small 
(e.g. <2 Ha) green spaces along stream corridors. 

We assessed variation in host-seeking abundance of mos- 
quitoes along a vertical gradient at four trapping locations 
in Foggy Bottom, DC between July and August 2010. We 
hung single trap-lines consisting of three CDC light traps 
baited with C0 2 at varying heights: 1.5, 6, and 10-15 me- 
ters with the highest trap height dependent on tree height. 

Roosting height and mosquito feeding preferences 

We examined how roosting height might influence mos- 
quito feeding preferences by using mosquito bloodmeal 
analysis data from a previous field study [28] conducted 
at the same field sites, and compared those data with the 
roosting height for the species in our study. In the previ- 
ous study, the preference for a species was estimated as 
the fraction of blood meals from a species divided by its 
relative abundance, estimated using four to six unlimited 
radius point counts at each site conducted monthly [28]. 
We limited our analysis to the average roosting height 
of juvenile birds, which appeared to differ in roosting 
height from adults although sample sizes prevented a 
rigorous statistical analysis. Thus, our analysis was re- 
stricted to five (American robin, Northern cardinal, Gray 
catbird, Mourning dove, and House sparrow) of the 
seven species in our roosting study. We had insufficient 
data to include European starlings in this analysis, and 
Tufted titmice were not present at the two sites where 
we had sufficient numbers of blooded mosquitoes to es- 
timate mosquito feeding preferences. 
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Simulated roost experiment 

We attempted to test the hypothesis that avian hosts roost- 
ing near a mosquito trap decrease the number of host- 
seeking mosquitoes caught in a trap. We compared mos- 
quitoes caught in CDC light traps baited with a single 0.66 
Gallon Igloo cooler full of dry ice with traps surrounded by 
four additional dry ice filled coolers spaced evenly around 
the trap in a spiral staircase pattern at a distance of 0.5 - 
1.0 m from the light trap with one in each cardinal direc- 
tion. We conducted this experiment at Foggy Bottom, DC 
and Takoma Park, MD. The experimental design was 
chosen to mimic the host densities observed around 
mosquito traps at actual communal roost sites. We es- 
timated that the four additional coolers produced the 
same amount of C0 2 per hour as 372.4 ± 24.8 add- 
itional hosts, based on estimated field metabolic rate of 
American robins [54]. We used field metabolic rate ra- 
ther than basal metabolic rate because many birds at 
communal roosts were moderately active during the 
night, including vocalizing and flying between perches. 
We ran traps for one night as described and then alter- 
nated the location of the additional coolers of dry ice. 
We replicated this experiment across four trapping lo- 
cations and eight nights for a total of 32 trap nights. 

Correcting mosquito abundance estimates at communal 
roosts 

We attempted to correct for the influence nearby hosts had 
on mosquito trapping efficacy by increasing mosquito 
abundance estimates at communal roost sites based on the 
results from the artificial roosting experiment. We calcu- 
lated the ratio of roost sizes of communal roosts and the 
number of American robins that would give off the same 
amount of C0 2 as our artificial roost experiment. We then 
increased mosquito abundance estimates using this correc- 
tion to calculate the amount of mosquitoes we would 
expect to catch at a communal roost site if nearby birds did 
not affect trapping success. 

Statistical analyses 

We compared differences in host-seeking mosquito 
abundance at communal and non-communal roost 
sites, and variation in avian host roosting behavior 
using generalized mixed models formulated with the 
glmer and Imer functions in the lme4 package in R (v 
2.15). For the mosquito abundance estimates we 
treated each individual mosquito trap and trapping 
period as random effects, and roost type (whether the 
collection site was used as an avian communal roost) 
as a fixed effect. In analyses of roosting height we used 
the identity of each individual bird as a random effect 
and modeled species as a fixed effect. Similarly, in ana- 
lyzing variation in roost aggregation size among avian 
species the identity of individual birds was treated as a 



random effect and species and the interaction terms of 
species by Julian date and species by the square root of 
Julian date (to account for non-linearity) were included 
as fixed effects. We selected the best fitting mixed 
models using Akaike's information criterion (AIC), and 
determined the significance of model terms by likeli- 
hood ratio tests on models with and without that 
predictor. 

Results 

Avian roosting and mosquito host-seeking behavior 

We deployed 213 radio transmitters nearly four-fold more 
than previous studies on WNV host behavior [40,41,55], 
and quantified roosting behavior (roost height and group 
size) of 91 individuals of six species from July-September. 
Roosting height did not differ among sites (^ = 0.602, 
df =2, P = 0.74) but differed among species (Figure la; 
mixed-effects model with individual bird as a random 
effect, X* = 36.59, df = 6, P < 0.001). Tufted titmice roosted 
highest in the canopy (data was primarily from at a single 
site near Rockville, MD), followed by American robins, and 
house sparrows roosted at the lowest height. 

Cx. pipiens host-seeking mosquito abundance also 
increased with trap height (Figure lb, generalized linear 
mixed-effects model with Poisson distribution and log-link, 
X* = 24.51, df = 2, P < 0.001). Finally, feeding preferences of Cx. 
pipiens mosquitoes for each host increased with average roost- 
ing height at the same sites (Figure lc; F 1>4 = 22.21, P = 0.018). 

Roost size, in numbers of birds, differed among species 
(Tables 1 and 2; species effect J 2 = 107.89, df = 6, P < 0.001) 
and roost sizes of three of the seven species, American 
robins, European starlings, and house sparrows, increased 
from July to September (Figures 2a-c, Tables 1 and 2). 
In late fall, American robins occurred in the largest 
roosts, followed by European starlings and house spar- 
rows. Mourning doves also roosted in groups, although 
aggregations for this species were much smaller than 
other communal roosting species (Figure 2d). Gray 
catbirds, northern cardinals, and tufted titmice often 
roosted individually or in pairs (Figure 2d). 

Raw estimates of host-seeking mosquito abundance 
varied among two communal and two non- communal 
roost sites (generalized linear mixed-effects model with 
log-link; site effect: jf 2 = 9.85, df=3, P = 0.02) but, sur- 
prisingly, traps at American robin communal roost sites 
caught significantly fewer mosquitoes than traps at non- 
communal roosting areas (Figure 3a; generalized linear 
mixed-effects model with log-link, x* = 4.03, df = 1, P = 
0.04). The vector :host ratios at communal roost sites 
using these raw data on host seeking mosquito abun- 
dance were 0.0011 ± 0.0004 mosquitoes per American 
robin per night at the two communal roost sites and 
2.62 ± 0.52 mosquitoes per American robin per trap 
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Figure 1 Host roosting height, mosquito foraging height, and host-vector associations, a. The roosting height of seven bird species in 
Maryland and Washington DC Bold crossbars indicate the mean, boxes encompass ± 1 SE of the mean, and whisker lines represent 95% confidence 
bounds around the mean. Note that all but one E. starlings studied were adult birds whereas all other species studied were composed of juvenile birds 
only. b. Mean number of host-seeking Go pipiens mosquitoes collected per trap night (N = 48) at three traps along a single trap line at four trapping 
locations (separated by vertical bars) in Foggy Bottom, DC. c. Mosquito feeding preferences from an earlier study [28] plotted against average roost 
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Table 1 Results of fitting a generalized linear mixed-effects 
model with a Poisson distribution and a log-link to 
analyze patterns of avian roost size over time for 
seven species of bird 



Predictor 


Coefficient 


SE 


Z 


P 


Intercept (A. Robin) 


-1075.1 


15.9 


-67.52 


<0.001 


T. titmouse 


1454.6 


602.8 


2.41 


0.016 


E. starling 


397.9 


150.2 


2.65 


0.008 


G. catbird 


866.5 


238.4 


3.64 


0.000 


H. sparrow 


566.7 


171.3 


3.31 


0.001 


M. dove 


1 888.5 


301.2 


6.27 


<0.001 


N. cardinal 


978.8 


134.1 


7.3 


<0.001 


Julian date (A. Robin) 


-4.4 


0.1 


-66.06 


<0.001 


T. titmouseJulian date 


6.0 


2.6 


2.32 


0.021 


E. starling:Julian date 


1.7 


0.6 


2.72 


0.007 


G. catbirdJulian date 


3.5 


1.0 


3.52 


<0.001 


H. sparrowJulian date 


2.3 


0.7 


3.11 


0.002 


M. doveJulian date 


7.9 


1.3 


6.06 


<0.001 


N. cardinahJulian date 


4.0 


0.6 


7.07 


<0.001 


Julian date 0 5 (A. robin) 


137.5 


2.1 


66.97 


<0.001 


T. titmouseJulian date 0 5 


-186.7 


78.8 


-2.37 


0.018 


E. starling:Julian date 0 5 


-51.5 


19.1 


-2.69 


0.007 


G. catbirdJulian date 0 5 


-110.2 


30.7 


-3.59 


<0.001 


H. sparrowJulian date 05 


-72.2 


22.4 


-3.22 


0.001 


M. doveJulian date 0 5 


-244.3 


39.6 


-6.17 


<0.001 


N. cardinahJulian date 0 5 


-125.1 


17.4 


-7.2 


<0.001 



American robin was the reference level against which other effects 
were compared. 



night at non-communal roost sites, or more than two- 
thousand fold higher (Figure 3b). 

Trapping bias - effect of hosts near traps 

Placing additional coolers of C0 2 near mosquito traps to 
simulate birds roosting nearby reduced the number of 
mosquitoes caught by 57% ±12% (Figure 4; generalized 
linear mixed-effects model with a Poisson distribution 
and a log-link, / = 298.04, df = 1, P< 0.001). We used 
this and the approximate number of birds that would 
produce the same amount of C0 2 (372.4 ± 24.8; see 
Methods), to adjust for the effect of additional hosts 
near mosquito traps at communal roost sites. This ad- 
justment increased mosquito abundance estimates from 
1.23 ± 0.4 to 3.6 ± 1.1 mosquitoes per trap night at com- 
munal roosts. With these adjusted estimates, there was 
no longer a significant difference in mosquito abundance 
between communal roosts and non-communal roosting 
areas (Figure 4; generalized linear mixed-effects model 
with a Poisson distribution and a log-link, j 2 = 0.0372, 
df=l, P = 0.85). However, the estimated vectonhost 



Table 2 Comparison of models to explain variation 
among species and over time in roost size using Akaike's 
Information Criterion (AIC), AAIC, Akaike weights [(o,) 



Model 


AIC 


AAIC 


61/ 


K 


Species*Jdate + Species*Jdate a5 


16694 


0 


1 


12 


Species + Jdate 


16965 


271 


0 


9 


Jdate + Jdate 05 


17135 


441 


0 


3 


Species*Jdate 


22144 


5450 


0 


10 


Species + Jdate + Jdate 05 


22266 


5572 


0 


10 


Jdate 


22418 


5724 


0 


2 


Species 


41269 


24575 


0 


8 


Intercept 


41373 


24679 


0 


1 



^denotes an interaction between two model parameters. 
Parameters in the model include linear and nonlinear time components 
(Jdate is the Julian date or day of the year, with January 1 = 1). K is the 
number of parameters in the model. 



ratios using adjusted mosquito abundance at roosts were 
still more than 850 times lower at communal roost sites 
(0.003 ± 0.001 mosquitoes per American robin) than at sites 
not used for communal roosting (2.62 ±0.52 mosquitoes 
per American robin, Figure 3b). 

Discussion 

Our results suggest two general aspects of host behavior, 
sociality and fine-scale habitat use, which have received 
relatively little research effort [40,41,55], influence pat- 
terns of vector-host interactions, and could thereby alter 
both transmission of and exposure to vector-borne path- 
ogens. Three species, American robins, European star- 
lings, and house sparrows, all formed communal roosts 
as mosquito abundance increased and infection preva- 
lence with WNV in this region and others subsequently 
decreased [29,30,40,41,56]. However, although large ag- 
gregations of birds produce substantial amounts of C0 2 
that might be expected to attract many host-seeking 
mosquitoes, per capita mosquito abundance at roost 
sites was two orders of magnitude lower than at non- 
communal roost sites. Even though we were not able to 
fully adjust our estimates of mosquito abundance for the 
distraction effect of birds roosting near mosquito traps 
(additional coolers of dry ice did not capture or satiate 
foraging mosquitoes, and did not include other odors 
emitted by birds), the enormous reduction in vector :host 
ratio suggests that sociality in these communally roost- 
ing species substantially reduces their contact with biting 
mosquito vectors. 

Reduction in average contact rates with biting vectors 
would reduce exposure to vector-borne pathogens. This 
would represent an additional fitness benefit conferred 
from communal roosting that complements other better 
studied effects such as protection from predators and 
centers for sharing foraging information [57]. Commu- 
nal roosting behavior in American robins, European 
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starlings and house sparrows occurred well before the 
arrival of WNV to North America and it is unclear how 
important the influence of other arboviruses such as St. 
Louis encephalitis virus has been in the evolution of 
communal roosting behavior. However, the introduction 
of WNV could have altered the selective pressure and 



timing of selection on communal roosting. The density 
of WNV-infected mosquitoes in the study area peaks in 
late July - early August [29], when communal roosts 
begin forming at our sites and elsewhere [40,41]. Given 
that approximately 12% of American robins infected 
with WNV die from infection [58], the selective pressure 
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Figure 3 Host-seeking mosquito abundance at communal and 
non-communal roost sites, a. Mean number of host-seeking Cx. 
pipiens collected per trap night from August through September 
2010 across (A) two non-communal roost sites and (B) two American 
robin communal roosts nearTakoma Park, MD. Adjusted mosquito 
abundance estimates for communal roost sites (C) account for the 
effect hosts have on mosquito trapping efficacy, b. Per capita number of 
host seeking Culex mosquitoes per night at (A) two non-communal roost 
sites and (B) two American robin communal roosts near Takoma Park, 
MD. Adjusted mosquito:host ratio estimates for communal roost sites (C) 
account for the effect hosts have on mosquito trapping efficacy. Error 
bars encompass ± 1 SE of the mean. 



of WNV on American robins to form communal roosts 
could be non-trivial [59] . 

Our results suggest that habitat use, - specifically for- 
aging and roosting height, by vectors and hosts, respect- 
ively, also plays a significant role in host-vector contact. 



As in previous studies [36-39], we found that host- 
seeking Culex abundance was nearly four times higher 
at mid-level traps and six times higher at high-level traps 
compared to traps placed 1.5 m from the ground. This 
suggests that mosquito-biting rates would be lowest on 
host species roosting lower in the canopy, and we found 
preliminary evidence of a correlation between mosquito 
feeding preference and avian roosting height. Higher 
host-seeking mosquito abundance in the canopy could 
result from C0 2 plumes from both traps and birds in 
the canopy spreading further than plumes from birds or 
traps at lower heights. House sparrows, which roost 
quite low (often in shrubs) may experience reduced 
mosquito contact rates due to mosquito foraging heights 
which could partially provide an explanation for the very 
low feeding preferences or host utilization on this spe- 
cies in many studies [27,28,30]. Further work on add- 
itional species and in other regions will determine 
whether roosting height is a consistent predictor of 
mosquito feeding preference. 

Previous work has shown that feeding of Cx. pipiens 
on a preferred avian host, American robins, decreases in 
the fall [27,28,30,44] and this decrease has been associ- 
ated with an increase in human feeding and human 
cases of WNV in one region [29]. This feeding shift was 
hypothesized to result from the late-summer dispersal of 
American robins following breeding. Our results suggest 
that the movement of American robins to communal 
roosts during this period might make them less available 
to host-seeking mosquitoes in urban and residential 
areas away from roosts and account for the reduced 
feeding on this species. The period when birds start 
forming communal roosts also coincides with a decrease 
in WNV infection prevalence in mosquitoes [29]. A de- 
crease in transmission at this time would be expected 
based our results that communal roosting reduces 
vector-host ratios on the most important species in 
WNV amplification in this and other regions, American 
robins [26-30,44]. 

Our findings raise questions about previous assertions 
that large aggregations of American robins increase WNV 
transmission [40]. While Diuk-Wasser et al. found higher 
WNV prevalence in mosquitoes at communal American 
robin roost sites compared to non-roost sites, they found 
no difference in total Culex mosquito abundance among 
roost and non-roost sites, and the ratio of per capita abun- 
dance of Culex mosquitoes at roost sites and non-roost 
sites was similar to what we found (<800 fold lower at the 
communal roost). This raises the question of how commu- 
nal roosting could have increased transmission at roost sites 
if per capita mosquito abundance was much lower. Two 
possibilities are: 1) that other non- roost related factors dif- 
fered among sites (roost sites are not randomly selected by 
birds or scientists) and 2) biting at roosts could be focused 
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on a very small subset of the birds at the communal roost 
resulting in a higher effective vectonhost ratio than one 
considering all birds in the roost. Additional work examin- 
ing WNV transmission ecology at a larger number of com- 
munal roost sites compared to sites not used for communal 
roosting would help answer this question, especially with 
mosquito trapping and testing for WNV occurring before 
roost formation. These studies should also account for the 
biases we have identified from mosquitoes being differen- 
tially attracted to traps at communal roost sites. 

The distribution of mosquitoes attempting to feed on 
individual birds within a roost also likely influences host 
defensive behavior and feeding success [32-34]. As a re- 
sult, communal roosting by birds likely disrupts the opti- 
mal distribution of mosquitoes feeding on each bird, 
because too many mosquitoes per bird are feeding at 
non-communal roost sites, likely leading to increased 
host defensive behavior and reduced mosquito feeding 
success [14]. 

Conclusions 

In summary, our results demonstrate how host aggrega- 
tion and habitat use and vector searching behavior influ- 
ence contacts between vectors and hosts. Differential 
habitat use by different host species and spatial hetero- 
geneity in vector host-seeking behavior partly explains 



why vectors rarely feed on hosts in proportion to their 
abundance. Our results also demonstrate the importance 
of quantifying the per capita rate of contact between 
vectors and communally roosting hosts and underscore 
the need to consider the spatial and among-individual 
heterogeneities in the disease ecology landscape. 
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